Optical bistability involving planar metamaterial with broken structural symmetry 
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We report on a bistable light transmission through a planar metamaterial composed of a metal 
pattern of weakly asymmetric elements placed on a nonlinear substrate. Such structure bears the 
Fano-like sharp resonance response of a trapped-mode excitation. The feedback required for bista- 
bility is provided by the coupling between the strong antiphased trapped-mode-resonance currents 
excited on the metal elements and the intensity of inner field in the nonlinear substrate. 
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The optical transmission bistability is a phenomenon 
whereby a system changes its transmission from one value 
to another in response to properties of light passing 
through The switch between the two stable states 
of transmission is usually induced by the intensity of the 
incident light, and, luckily, the switch may exhibit a hys- 
teresis. The light intensity driven systems must have an 
intrinsic nonlinear susceptibility and feedback to provide 
optical bistability. This effect can be used to realize all- 
optical switches, limiters and logic gates. The issue of 
the day is to reduce the size of such all-optical devices, 
decrease their switching times and the required intensity 
of light. 

Typically, all-optical devices include a Fabry- Perot res- 
onator to provide feedback 0; although the systems 
without such resonator can be used too. An example 
of such structures is the nonlinear photonic crystal mi- 
croresonators |3j, |4| and the devices based on the surface 
plasmon-polaritons in metal nanostructures [5|, |6| . They 
offer a unique mechanism for confining light, giving rise 
to a combination of high quality factors and small modal 
volumes that is helpful to stronger nonlinear effects and 
reduce the size of the all-optical devices. 

There is another promising way to develop the small- 
size all-optical switching devices without electromagnetic 
cavities and surface plasmons field enhancement. Here 
the case in point is the planar metamaterials (also known 
as metafilms) with active constituents. Typically, these 
systems are surface structures which consist of some 
metal or dielectric resonance elements arranged as a peri- 
odic array and placed on a layer thin in comparison to the 
wavelength. One way to obtain the nonlinear response 
of a planar metamaterial is to introduce some nonlin- 
ear individual resonance elements in the structure (see 
for instance (tJ)- Thus in (gj, the elements are made 
nonlinear and tunable via the insertion of diodes with a 
voltage-controlled capacitance. However, in the optical 
range the manufacturing of such structures is associated 
with considerable technological difficulties. Another sim- 
pler way is to arrange the proper resonance elements on 
a nonlinear substrate. 

The main feature of the planar metamaterials, essential 



for the optical switching applications, is a resonance char- 
acter of their transmission and reflection spectra. The 
usual resonance field enhancement inside a planar meta- 
material may be extremely enlarged by involving struc- 
tures which bear so-called trapped or dark modes 0, [l(| ■ 
The excitation of high-quality-factor trapped mode res- 
onances in planar double-periodic structures with a bro- 
ken symmetry was shown both theoretically llj, [l2[ and 
experimentally 13] in microwaves. In particular, these 
typical peak-and-trough Fano spectral profile resonances 
are excited in the periodic structure consisted of asym- 
metrically split metal rings. A small asymmetry of the 
metal elements of such structure results in the excita- 
tion of the strong mode of anti-phased currents, which 
provides low radiation losses and therefore high Q-factor 
resonances. Recently, the trapped mode resonances were 
investigated in similar planar structures in the near-IR 
range [14J. It was shown that the special choice of ge- 
ometry parameters of the structure enables to increase 
the Q-factor of the trapped mode resonance by several 
times in comparison with the ordinary plasmon-polariton 
resonance. Such high-Q resonance regime is promised to 
observe a bistability in the near-IR range, if the structure 
is to include a nonlinear material. 

In this report, we propose and study an all-optical 
switching device based on the bistability in a planar 
metamaterial made of complex shaped resonance parti- 
cles placed on a substrate of nonlinear material in the 
regime of a trapped mode excitation. 

The studied structure consists of the identical subwave- 
length metal inclusions in the form of asymmetrically 
split rings (ASRs) arranged in a periodic array and placed 
on a thin nonlinear dielectric substrate (see Fig.[lJ. Each 
ASR contains two identical strip elements opposite to one 
another. The right-hand split between the strips ip± is 
a little different from the left-hand one <p2, so that the 
square unit cell is asymmetric with regard to the y-axis. 
A normally incident plane wave polarized transversely to 
the array symmetry axis (x-axis) can excite the trapped 
mode resonances. Therefore we consider the structure 
excitation by y-polarized plane wave. Suppose, that the 
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ing A = 1 v/m the magnitude of the current I averaged 
along the single element can be determined as a function 
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FIG. 1: (Color online) Fragment of the planar metamaterial 
and its elementary cell. The size of the square translation cell 
is d = d x — d y = 900 nm. The asymmetrically split metal 
rings are placed on the top face of h = O.fd thick nonlinear 
dielectric layer. The radius and the width of metal rings are 
a = OAd and 2w = 0.06d respectively. The angle sizes of ring 
splits are <pi = 15° and if 2 = 25°. 



incident held has an amplitude A and a frequency u>. 

As usual, actual optical array structures are mounted 
into a dielectric host or placed on a thick dielectric sub- 
strate. Two things happen with the trapped mode reso- 
nance frequency loq , when a dielectric substrate is added 
next to the structure. First, the resonant frequency ujq 
changes. If the structure is completely inside the dielec- 
tric material with relative dielectric constant e r , the res- 
onant frequency will reduce in frequency with the fac- 
tor yfe^ [I5j]. If the structure is placed between dielec- 
tric slabs of finite thickness, the resonant frequency will 
change to somewhere between ojq and wo/y^r- And, at 
last, if we have a dielectric only to one side of the struc- 
ture, the resonant frequency shifts downward to be near 
luq/ yj (e r + l)/2. Second, the finite thickness of a practi- 
cal substrate (of about 0.3-0.5 mm) results in the appear- 
ance of some additional interference resonances. How- 
ever, as shown in the trapped mode resonance can 
be easily separated against the interference resonances 
background. Thus the presence of the thick dielectric 
substrate has no effect on the trapped mode resonance 
qualitatively; and we consider further a "free standing" 
structure configuration of a metal array placed on a thin 
nonlinear substrate, in order to simplify explanation. 

The algorithm based on the method of moments was 
proposed earlier [H| to study the resonance nature of 
the structure response, under the assumption of such a 
small amplitude A that the dependence of the substrate 
permittivity e on the field intensity is infinitesimal. 

The algorithm requires that, at the first step, the sur- 
face current induced in the metal pattern by the field of 
the incident wave is to be calculated. The metal pattern 
is treated as a perfect conductor, while the substrate is 
assumed to be a lossy dielectric. In particular, assum- 



I = Q(u,e) 



(1) 



At the second step, the found surface current distribu- 
tion is used to calculate the transmission and reflection 
coefficients as t — t(w,e), r = r(uj 7 e). 

To introduce the nonlinearity (the third-order Kerr- 
effect), let us assume that the permittivity s of the sub- 
strate depends on the intensity of the electromagnetic 
field inside it. In our approximate approach to the nonlin- 
ear problem solution, first assume that the inner intensity 
is directly proportional to the square of the current mag- 
nitude averaged over a metal pattern extent. Secondly, in 
view of the smallness of the translation cell of the array, 
we suppose that the nonlinear substrate remains to be a 
homogeneous dielectric slab under intensive light. Thus, 
the intensity-dependent permittivity of the substrate is 
given further as 

e = e 1 +e 2 \I\ 2 . 

Note, since the substrate permittivity is proportional to 
the current value in the metal pattern, the nonlinear- 
ity effect reaches its maximum under the ASR resonance 
condition. 

If the amplitude A of the incident field differs from the 
unity, the appropriate average current magnitude for a 
given e can be found using (JT|) as 



I = A-Q(u,e). 



(2) 



Since the substrate permittivity e depends on the av- 
erage current value /, the relation ([2]) can be rewritten 
as follows 



1 = A-Q(u, £l +e 2 \I\ 2 ). 



(3) 



The expression (|3]) is a nonlinear equation related to 
the average current value in the metal pattern. The inci- 
dent field magnitude is a parameter of the equation ([3]) . 
At a fixed frequency ui, the solution of this equation is 
the average current value dependent on the magnitude of 
the incident field I = 1(A), where the function 1(A) is 
presumably multivalued. 

On the basis of the current 1(A) found by a numer- 
ical solution of the equation @, it is possible now to 
determine the permittivity of the nonlinear substrate 
e = E\ + £2 1 1(A) 1 2 and to calculate the reflection and 
transmission coefficients 

* = t(a J ,e 1 +e 2 \I(A)\ 2 ), r = r(«, £1 + e 2 \I(A)\ 2 ), 

as the functions of the magnitude of the incident field. 

At first we consider the transmission through the a rray 
of ASRs placed on a linear substrate (see Fig. [5]) plj - flij ]. 
If a normal incident wave is polarized in y-direction, at 
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FIG. 2: (Color online) The frequency dependences of the mag- 
nitude of the transmission coefficient (a) and the square of 
the current magnitude (a.u.) averaged along the metal pat- 
tern (b) in the case of the linear permittivity (52 = 0) of the 
substrate. 



FIG. 3: (Color online) The square of the current magnitude 
in a.u. averaged along split ring (a) and the magnitude of the 
transmission coefficient (b) versus the incident field magni- 
tude in the case of the nonlinear permittivity (ei = 4 + 0.02i, 
£2 = 5 x 1CP 3 ) of the substrate. The arrows indicate the 
bistable loop. 



the dimensionless frequency nearly ae = d/X ~ 0.3, a 
sharp reflection resonance occurs. This resonance corre- 
sponds to the excitation of a trapped mode because equal 
and opposite directed currents in the two arcs of each 
complex particle of array radiate a little in free space. 
The resonance has a high quality factor, and the current 
magnitude reaches the maximum at this frequency. As 
the permittivity of the substrate increases, the resonance 
frequency shifts to low values. Note, if the incident field 
is x-polarized or the splits between the strips are the 
same, only a symmetric current mode is excited. The 
corresponding resonance has a low quality factor, and 
this configuration is not suitable to observe a bistability 
because the current magnitude is small at this frequency. 

Suppose that the trapped mode resonance frequency 
is slightly higher than the incident field frequency. As 
the intensity of the incident field rises, the magnitude of 
currents on the metal elements increases. This leads to 
increasing the field strength inside the substrate and its 
permittivity as well. As a result, the frequency of the 
resonant mode decreases and shifts toward the frequency 



of incident wave, which, in turn, enhances further the 
coupling between the current modes and the inner field 
intensity in the nonlinear substrate. This positive feed- 
back increases the slope of the rising edge of the trans- 
mission spectrum, as compared to the linear case. As the 
frequency extends beyond the resonant mode frequency, 
the inner field magnitude in the substrate decreases and 
the permittivity goes back towards its linear level, and 
this negative feedback keeps the resonant frequency close 
to the incident field frequency. 

The curves in Fig. [3^ are a plot of the square of the 
current magnitude averaged along ASR element versus 
the incident field magnitude at fixed frequencies. As 
an example, let us consider the incident field frequency 
a; = 0.301 (solid line). As the incident field amplitude 
increases, the current magnitude (and proportionally the 
inner field intensity) gradually increases along the bot- 
tom branch of the curve, until it reaches about \I\ 2 f=s 2. 
At this point, the current magnitude jumps to around 
\I\ 2 « 5.5 due to the instability of the system at the inte- 
rior branch of the curve. This transition is shown in the 
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FIG. 4: (Color online) Frequency dependences of the square 
of the current magnitude averaged along split ring (on the 
logarithmic scale) in a.u. (a) and the magnitude of the trans- 
mission coefficient (b) in the case of the nonlinear permittivity 
(ei = 4 + 0.02i, e 2 = 5 x 10~ 3 ) of the substrate. 



figure with an arrow directed upward. The incident field 
magnitude decreasing results in flowing the square of the 
current magnitude along the upper branch of the curve 
down to \I\ 2 ps 4.5, where it drops to a value of about 
\I\ 2 fa 0.5. Similarly, this transition is displayed as an ar- 
row directed downward. The dramatic current variation 
from high to low level produces a switching from high to 
low transmission (see Fig.[3jD). 

The frequency dependence of the transmission coeffi- 
cient magnitude manifests also some impressive discon- 
tinuous switches to different values of transmission, as 
the frequency increases and decreases in the resonance 
range for the sufficiently large intensity of the incident 
wave. The shifting of the peak of the resonance and the 
onset of a bistable transmittivity through the ASR struc- 
ture is similar to that of the reflection from a Fabry- Perot 
cavity (Fig. UK). However, the trapped- mode resonance 
is Fano-shaped rather than the Lorentzian, as is the 
characteristic of ID Fabry-Perot cavities. This Fano res- 
onance can lead to a peculiar transmission spectra and 
bistable behavior [13]. In particular, the transmission 
resonance of the ASR structure may loop back on them- 



selves (Fig. [4)3). Such a form of the hysteresis loop of 
transmittivity is similar to the behavior of the reflection 
spectra of a photonic crystal in the presence of a non- 
resonant downstream scattering source in microcavities 
17]; where the Fano resonances arise due to the inter- 
ference between two (or more) different scattering path- 
ways. From the analogy between these two systems, it 
can be concluded that the peculiar bistable loop in the 
transmission is the result of a superposition of resonant 
and nonresonant contributions. 

In conclusion, a planar nonlinear metamaterial com- 
posed of a metal array placed on a nonlinear dielectric 
substrate, which bears a sharp resonance response by an 
excitation of a trapped mode due to a broken symme- 
try of the pattern is a challenging object for all-optical 
switching applications. 

This work was completed with the support of the Na- 
tional Academy of Sciences of Ukraine by the Program 
" Nanotechnologies and Nanomaterials" , the Project 
no. 1.1.3.17. 
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